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In the two-dimensional case the transition between ordinary and topological insulator states can 
be described by a massive Dirac model with the mass term changing its sign at the transition point. 
We theoretically investigate how such a transition manifests itself in resonant transport via localized 
helical edge states. The resonance occurs in the middle of the band gap due to a zero edge-state 
mode which is protected by the time-reversal symmetry, also when coupled to the conducting leads. 
We obtain the explicit dependence of the resonant conductance on the mass parameter and an 
external magnetic field. The proposal may be of practical use, allowing one to determine the orbital 
g-factor of helical edge states in two-dimensional topological insulators. 

PACS numbers: 



I. INTRODUCTION 

Two- (2D) and three-dimensional (3D) topological in- 
sulators (TIs) are in the focus of current research (see, 
recent reviews^— and references therein). They differ 
from ordinary band insulators by the presence of pro- 
tected edge (in 2D) or surface (in 3D) states that sup- 
port electric current and, for that reason, may have some 
application potential. The 2D TIs are of special impor- 
tance because they exhibit the quantum spin Hall (QSH) 
stated— that, unlike the quantum Hall (QH) states, does 
not require any magnetic field and is characterized by 
time-reversal invariant gapless edge modes, while the 
bulk states are fully gapped. Such edge modes, called 
sometimes "helical", consist of a pair of channels that 
propagate in the opposite directions for opposite spin di- 
rections on the same edge. Another interesting feature of 
the 2D TIs, which distinguishes them from 3D TIs, is that 
their low-energy physics can be described by a massive 
Dirac model££ In particular, the transition between the 
ordinary and topological insulating regimes corresponds 
to inversion of the effective relativistic mass of the Dirac- 
like fermions. At the transition the system behaves as a 
gas of massless Dirac fermions, 8 

The first 2D TI was realized experimentally in HgTe 
quantum wellsji&i In these experiments, the QSH regime 
was detected by measuring the longitudinal conductance 
of two spin channels propagating in the same direction 
on opposite edges of the sample. This finding was fur- 
ther substantiated by the observed suppression of the 
edge transport in a magnetic fieldji^ which breaks the 
time-reversal symmetry of the QSH state, thus revealing 
the helical nature of edge channels. Theoretically, the 
transport in QSH insulators has been analyzed in both 
disordered^— and ballistic^ - — regimes. 

In the present paper, we discuss another transport 
regime - resonant tunneling - as an alternative means 
of probing helical edge states in 2D TIs. We consider 
edge states localized on a finite-length boundary coupled 
by tunneling to two metallic leads (see also Fig. [IJ. The 



operation of such a device relies on the fact that gap- 
less helical edge states are protected by the time-reversal 
symmetry!^ Therefore, for any tunneling coupling that 
preserves the time-reversal invariance the edge spectrum 
contains a doubly degenerate zero mode in the middle 
of the band gap, which provides the resonant level for 
tunneling between the leads. If, upon the gap inversion, 
the systems goes into the ordinary insulating regime, the 
zero mode disappears and the transport through the de- 
vice is suppressed. This can be described by the following 
formula for the tunneling conductance: 



g{M,B) = - 



iRlL 



_{M + \M\ + fiBf + (2E±my 



iRlL 

(M + \M\ - [IB)* + (i*p±y 



(1) 



where M is the mass (gap) in the Dirac model 5 and jl/^ 
and Jr/H are the tunneling rates between the edge and 
the left (L) and right (R) contacts. The two terms in 
Eq. (JXJ) arise from lifting the Kramers degeneracy of the 
zero mode in an external magnetic field B, which involves 
the effective orbital magnetic moment of the edge states 
\i = \e\hv 2 /2c\M\ (v is their velocity). 

The transition between the topological and ordinary 
insulator states manifests itself in the strong dependence 
of conductance JT]) on the sign of M (see, also, Fig. QJl). 
In the topologically nontrivial state with M < 0, the 
B-ficld dependence of conductance ([TJ has a resonant 
(Lorenzian) line-shape around B = 0, while in the or- 
dinary insulator regime with M > the conductance is 
small and field-independent, g « (e 2 / h)"f l / 2M 2 for 
1l,r < \M\ and [iB < 2\M\. 

Apart from the principal possibility to detect the edge 
zero mode, the resonant magnetoconductance ([T]) can 
be used for determining the effective orbital g-factor 
of the helical edge states in a given material. Such a 
need exists, e.g., in magnetotransport studies of HgTe 
quantum wellsi For the typical velocity in HgTe wells, 



2 



v = 5.5 x 10 5 ms _1 fi we estimate \i (in units of Bohr 
magneton hb) as h/hb ~ 1700meV/|M|, where M is 
measured in meV. For the realizable values of M ~ 10 
meV, the orbital g-factor is quite large ~ 170. 

The subsequent sections describe the details of our the- 
oretical analysis. In Sec.|H]we solve the boundary prob- 
lem for helical edge states in a QSH insulator subject 
to a finite magnetic field. Section IIIII presents a more 
general formula for the resonant magnetoconductance in 
finite magnetic fields, from which we derive Eq. ([1]), and 
finally Sec. [IV] concludes the paper. 




insulator 



L 



. 1 




y 






Contact (1 




QSH 


insulator 

z 

®B 


Y* 







insulator 



I 



II. QUANTUM SPIN HALL INSULATOR IN A 
MAGNETIC FIELD 



A. Boundary problem 

We begin by analyzing the edge states on a boundary 
of a QSH insulator of finite length L, imposing symmetric 
boundary conditions at the ends x = ±L/2 (see, Fig.QJ,). 
The symmetric boundary conditions result in a vanishing 
particle current density at x = ±L/2, thus providing a 
confinement for Dirac fermions in a QSH insulator (see 
also Appendix |A"|) . For concreteness, we use the effective 
4-band model derived in Ref. d for HgTe quantum wells. 
In this approach one works in the basis of the four states 
near the T (p = 0) point of the Brillouin zone: |ei+), 
\hi+), |ei — ), and \h±— ), where e\ and h\ are the s-like 
electron and p-like hole QW subbands, respectively. The 
index r = ± accounts for the spin degree of freedom. 
The effective two-dimensional Hamiltonian can be ap- 
proximated by a diagonal matrix in t spaced 



H 



'p 








d p <T, dp 



{vp x ,-vp y ,M). (2) 



where Pauli matrices u x ,y,z act in subband space, v sa 
5.5 x 10 5 ms _1 is the effective velocity^, and M determines 
the bulk band gap at p = 0. In Eq. ([2]) we omit terms 
cx p 2 since our main results will hold for the vicinity of 
the p = point. Up to a unitary transformation, Eq. @ 
is equivalent to a massive Dirac Hamiltonian 



H D = vr z ap + Mt z g z , 



(3) 



t z is the Pauli matrix in spin space. We will work with 
the corresponding retarded Green's function defined by 



{eI-H D }G{v,v') = I8{T-v% 



(4) 



where 



p = -7wV-eA(r)/c, A(r) = {-By, 0, 0). (5) 

A(r) is the vector potential of an external magnetic field 
B, and / = r a = diag(l, 1, 1, 1). 

The system size in the y direction is considered much 
larger than the characteristic decay length of the edge 
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FIG. 1: (a) Counter-propagating (helical) edge states at one 
of the boundaries (y — 0) of a QSH insulator, (b) and (c) 
Schematics of resonant tunneling via localized edge states on 
finite-length linear and circular boundaries in a QSH insula- 
tor. In (c) the edge states appear on the outer boundary of 
an island (e.g. suitably chosen ordinary insulator or a hole) 
surrounded by a QSH insulator. 7i/ft and 7ij/A are tunnel- 
ing rates between the edge and the left (L) and right (R) 
contacts, (d) Magnetoconductance fU for topologically non- 
trivial (M < 0) and ordinary (M > 0) insulator states for 
7l = -jr = 0.02 \M\. Inset: Enlarged view of g{B) for M > 0. 



states, ~ hv/\M\, allowing us to focus on one of the edges 
{y = in Fig. QJi). At y — we impose the condition^ 



G(r, r')\ y=0 = T ®a x G(r, r')\ y=0 . 



(6) 



It can be obtained by introducing a large mass term 
{M — > +oo) outside the physical area of the system. 
Therefore, with M < in the bulk of the system the 
boundary condition ([6]) is equivalent to a mass domain 
wall which is responsible for the appearance of the edge 
states in a 2D QSH insulator in close analogy with the 
3D situation considered in Ref. One can check that 
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the boundary condition (J6j) indeed ensures vanishing of 
the normal component of the particle current: 

j y (x,Q) = ^(x,0)T 3 <S><7 y ij(x,0) (7) 
= i>\x, 0)t (g) <t x (t 3 <g) <t v )to <E> a x ip(x,0) 
= -V' t (x,0)r 3 <8> a y ip(x,Q) = -j y (x,0) = 0, 

where we switched to the creation and annihilation op- 
erators subject to boundary condition (|6]): i/j(x,0) = 
to <g> Ox ip(x, 0) and tp^(x, 0) = tp^(x, 0)tq ® a x . 



B. Green's function solution 

The solution of this boundary problem can be ob- 
tained, following the same steps as in the problem of 
the edge states in graphene . 20 ! 21 The Green's function is 
block-diagonal in Kramers partner (r) space: 



G 



G+ 
G_ 



G\\\ T Gi2\ 7 



G 2 l| r G 



22|t 



±, (8) 



where G T is a matrix in space of the two QW subbands 
(or, generally, in space of the upper and lower compo- 
nents of the Weyl spinor), labeled by indices 1 and 2. 
Writing Eq. (|4]) in components, it is easy to express the 
off-diagonal elements G 12 | r and G 21 | r in terms of the di- 
agonal ones as follows 



G T 



Gn| T 

: Lri 



vp- n 

Cz 22|r 



(Y - Y ) 2 



Q Y 2 



- 0,2 



e-Mr , 

G 22 | Tn = fe2 „.2 / \ "v* ~Y ), 



H 2 V 2 /X 



a 2 



M 2 - e 2 sgn(eB) 



(hv/Xf 



(11) 



where position Y = y/X is measured in units of magnetic 
length A = (ch/2\eB\) 1 ^ 2 , and F = -2A fc„ sgn(eB) is 
the center of the oscillator. The boundary conditions for 
Eqs. ([TUl) and (fTTj) are derived from Eq. ([6]) by writing it 
in components and using Eq. 



dG lllrn (Y,Y>) 



re + M y sgn(eB) 
7w/A + 2 



«iGn|rn(0,r'), (12) 
(13) 



5G 22 |Tn(^ Y ) 



ay 



k 2 



= K2 G 22 | T „(0,F'), (14) 
re - M y sgn(eS) 



fe/A 



(15) 



We seek the solution to Eq. (fTU|) in the form of the 
linear combination: 



1|t 



G 



22|7 



TE+M 

where for Gu\ T and G 22 | T we have 

[e 2 - M 2 - v 2 p-p + ]G nlr = (e + Mt)5(t - r'), 
[e 2 - M 2 - v 2 P+P _]G 22lT = {e- Mr)5{v - r>). 

Expanding 



G T (r,r') = ]T G Tn (y,y')e ik ^-^/L, k n =2m/L, 



,p±=p x ± ip y , (9) G n[Tn (Y, Y') = G^ llTn (Y, Y') + AxU{a x ,Y - Y ). (16) 



The first term is the Green's function of the unbounded 
system (source term), 



r , _ r , U( ai ,Y -Y )U(a u -Y' + Y )\y>y>, 

l*U\m ~ ^ 1 UfaY' - Y Q )U{a U -Y + Y Q )\y,>y, 



C X 



A(e + Mr)r( ai + 1/2) 
2^h 2 v 2 



(17) 



we obtain ordinary differential equations for Gn| rn (y, y') 
and G 22 | T „(y,?/') : 



d 2 

dY 2 



(Y-Y f 



a i 



Gn| Tr . 



e + Mr 

h 2 v 2 /x 



S(Y - Y') 



M 2 



a i 



sgn(eB) 



(fe/A) 2 



(10) 
I 



where U(ai,Y — Yq) is the parabolic cylinder function 
(see, e.g. Ref. \22h and T(ai + 1/2) is Euler's gamma 
function. The second term in Eq. (TT6l) is the solution of 
the corresponding homogeneous equation, decaying for 
Y — ^ oo . The coefficient A\ can be found from boundary 
condition (TT2"1) . which finally yields 



We have also used the recurrence relation U'(a, x) — (x/2)U(a, x) — U(a— 1, x)M- G 22 | T „(Y, Y') is obtained similarly 
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and has the form: 



G 22lrn (Y,Y>) = G%, n (Y,Y')-C. 



(Y /2 + K2 )U(a 2 , Y ) - U(a 2 - 1, Y Q ) 
! (Yq/2 + K 2 )U(a 2 ,Y ) + U(a 2 - 1, Y ) 



U{a 2 ,Y-Y )U(a 2 ,Y' -Y ), (19) 



with the source term 



LT 22|rn 



Co 



= c 2 
xu- 



U{a 2 , Y - Y )U(a 2 , —Y' + Y )\ y >y> 
U{a 2 , Y> - Y )U(a 2 , —Y + Y )\ Y ,> Y 

Mr )T(a 2 + 1/2) 



(20) 



C. Transition between ordinary and topological 
regimes in weak magnetic fields 

Let us now demonstrate how the Green's function for- 
malism describes the transition between ordinary and 
topological regimes. We will assume that energy e and 
magnetic field B are small in the sense that 



|e|<|M|, Hv/X<\M\. 



(21) 



This corresponds to the positive and large a\ [see, 
Eq. (|10p], when the parabolic cylinder functions in 
Eq. (JTSJ assume the form U(a, x) « 0F/(2 a / 2+1 / 4 r(3/4+ 
a/2))e~v* x Using then Stirling's expansion for T(x), 
we find 



Gii\ Tn (y,y') ' 
e + tM M 



+ 



\M\-*£n- 



2hv\M\M- 



\M\ -!fn- 



(22) 

-\M\(y+y')/hv 



fiB 



The nontrivial dependence of the edge (second) term on 
the mass sign describes the transition between the ordi- 
nary and topological insulator states. Indeed, the gapless 
edge modes appear only when M reverses its sign from 
positive to negative: 



hv 

— Tn + [irB, M < 0. 

1j 



(23) 



Unlike the conventional quantum Hall systems^ such 
edge states exist in a zero or weak magnetic field, when 
their localization length Kv/\M\ is determined by the 
band gap \M\ rather than the magnetic field £>. 24 The 
weak magnetic field results in a Zeeman-like term [itB 
in Eq. ([23]) where the effective magnetic moment /j, — 
\e\hv 2 /2c\M\ has an orbital origin: it is related to the 
magnetic flux through the finite width Hv/\M\ of the edge 
state. The magnetic field lifts the Kramers degeneracy 
e_ Tj _„ = e Ti „, splitting the Kramers partners by 2[iB in 
energy. 

In the general case, the analytical results for the edge- 
state spectrum are corroborated by the numerical solu- 
tion of equation (Y /2 + Ki)U(ai,Y ) + U(ai-l,Y ) = 0, 
which corresponds to the pole of Green's function (Q2 



For positive M we again find no subgap states, while for 
M < the solution is shown in Fig. [2] In a vanish- 
ing magnetic field (see, Fig. [5fx) , the spectrum contains 
a Kramers-degenerate zero mode with n = in the mid- 
dle of the band gap. The finite field lifts the Kramers 
degeneracy, opening a gap at n = (see, Fig. [^b). We 
emphasize that the gap opening is the consequence of the 
discreteness of the spectrum. For the continuum spec- 
trum (i.e. for L — > oo) the effect of the orbital magnetic 
field reduces to a shift of the crossing point of the r = ± 
branches. 17 



D. OD case 

Hereafter we focus on the OD case which admits a sim- 
pler analytical treatment at finite magnetic fields. This 
limit is realized when the level spacing hv/ L is much 
larger than both thermal activation energy fc^T and level 
broadening 7 (due to the coupling to metallic leads, see 
Sec. El: 



hv/L > ksT, 7. 



(24) 



Under such conditions, the transport at energy e = is 
determined by the presence (or absence) of the Kramers- 
degenerate zero mode with n — (see, Fig. [2]). Conse- 
quently, in Green's functions (fT5)l and p^|) one can set 
Y = and use U(a,0) = x A : /2 a / 2 + 1 / 4 r(3/4 + a/2) (see, 
Ref. [H), which yields 
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FIG. 2: Energies of spin-up e+, n and spin-down e_, n edge 
states (in meV) vs. wave number n for M = —5 meV, 
v = 5.5 x 10 5 ms _1 and L = 1 /im. (a) In a vanishing magnetic 
field the spectrum is gapless with the Kramers-degenerate 
zero mode at n — 0. (b) The degeneracy is lifted in a fi- 
nite field opening a gap at n — 0. For B — 60 mT the gap 
energy is approximately 1 meV. 
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Gn\r (Y,Y r ) ~ -C! 



A(e + Mt)B(1/2, 1/4 + ai/2) - sfimKvr 
A(e + Mr)S(l/2, 1/4 + oi/2) + ^Hvt 



U{ai,Y)U(ai,Y'), 



(25) 



G 22 | T o(r,r')~-^ 



A(e - Mr)B(l/2,l/4 + g 2 /2) + y/2nhvT 
A(e - Mt)B(1/2,1/4 + a 2 /2) - V^hvr 



U{a 2 ,Y)U{a 2l Y'), 



(26) 



where f?(l/2,l/4 + a/2) is Euler's beta function (we 
have omitted the bulk terms since they do not have the 
poles inside the band gap). Next, we integrate Eqs. (|2"5)) 
and (|26p over the remaining variable y and introduce 
the OD functions as Gn| T0 (e) = J °° Gh\ t q(Y, Y)dY and 



G 2 2|T-o(e) = Jo°° G 2 2\to(Y, Y)dY. The integrals can be 
expressed in terms of the digamma, ip, and gamma func- 
tions using the identity U 2 {Y)dY = vWf[VK 3 / 4 + 
a/2) - V>(l/4 + a/2)]/2r(l/2 + a) (see, Ref. M): 



_ A 2 (e + Mt)[i/;(3/4 + ai/2) - V>(l/4 + ai/2)] A(e + Afr)fl(l/2, 1/4 + ai /2) - V^Hvt 
11|r0l£j ~ 2^ 2 u 2 A(e + Mr)S(l/2,l/4 + ai/2) + V2wfer' 1 ' 



_ A 2 (e - Mt)[^(3/4 + ag/g) - V>(l/4 + a 2 /2)] A(e - Afr)fl(l/2, 1/4 + a 2 /2) + V27rfar 
22|r0l£j ^ 2nK 2 v 2 A(e-Afr)S(l/2, l/4 + o 2 /2)- V^hvr 



In the next section we will use these equations to calcu- 
late the resonant magnetoconductance through the zero 
edge mode. 
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FIG. 3: Conductance (J33J in units of 2e 2 /h versus gap M 
(meV) for B = 1 mT. Inset shows the same dependence for 
higher field B = 100 mT; v = 5.5 x 10 5 ms~ 1 and -y L = 7_r = 
0.3 meV. 



III. RESONANT MAGNETOTRANSPORT 

Transport through the localized helical edge states can 
be realized in the geometries shown in Figs. Qi> and c 2e 
where the edge is coupled to two metallic leads (L and R) 
via tunneling. Following the general approach of Ref. [27| 
we describe the coupling by the tunneling Hamiltonian, 



Hn 



E E 



d 



H.C.), 



(29) 



where operator c^ a (ck a ) creates (destroys) an electron 
with momentum k and spin a in lead L or R, and d\ lT 
and d nT are the creation and annihilation operators for 
a localized edge state with quantum numbers n and r. 
The tunneling matrix elements Vka.nr are assumed spin- 
independent. Further, we will work in the OD limit (|24|) 
where only the zero edge mode, described by d\ T and 
do T in Eq. (|29p , needs to be taken into account. Since 
the problem is similar to the tunneling through a reso- 
nant level, we can apply the known formula for the zero- 
temperature two-terminal conductance^ 



9 



e 2 ^IrIl 
h 7_r + 7l 



Trlm£?(e = 0), 



(30) 



where the tunneling energies Jl.r are expressed in terms 
of the tunneling matrix elements and Green's functions 
of the leads<21 In Eq. ([50]) Trim denotes the trace in 
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T®a space of the imaginary part of the retarded Green's 
function (/(e), which is related to the Green's func- 
tion G(e) calculated above through the Dyson equation, 
G(e) = G(e) + G(e)£0(e), involving the tunneling self- 
energy S. For metallic leads and in the absence of spin 
scattering, S is energy-independent and diagonal in T®a 
space: 



27 



-ryi", 7= ( 7b + 7l )/2. 



(31) 



The solution of the Dyson equation takes the form Q(e) = 
G(e + ij). This self-consistently accounts for the edge- 
spectrum broadening due to the presence of the leads. 
With the OD Green's functions (J^J) and ([!§, Eq. (j30)) 
reads 



!1 



or, explicitly, 



Im[G U | T0 (i7) + G aa , T0 (t7)],(32) 



g(M,B) = 



2e^ r r 1 (x)r 1 (x) 7fi7L 

h [(M + \M\T 1 (x)/ y ^) 2 + 7 2 ~ 
T 2 (x)T' 2 (x) 



(M+\M\T 2 (x)/^) 2 + 1 2 



= B i«/l s l 



(33) 



where B M = cM 2 / (4\e\hv 2 ) and the functions Ti :2 (x) 
and their derivatives T'^ 2 (x) are expressed in terms of 
the gamma functions: 

ri(a;)=r(a; + l/2)/r(a;), T 2 (x) = x/T^x). (34) 

Let us analyze the results for the magnetoconduc- 
tance. Equation (|33|) is valid for arbitrary ratio of B 
and the characteristic field, B M , related to the band gap 
M. Equation ([I]), presented in the introduction, fol- 
lows from Eq. ([33]) in the case of weak fields (or suffi- 
ciently large gap |M|) where x~ x — \B\/B M < 1. To see 
this we use Stirling's formula to obtain the expansions 
T^x) « - x^/S) and T 2 (x) w y/x{l + x~ l /%) and 

then insert them into Eq. (f3"3"|). 

As \M\ becomes smaller (i.e. x = B M /\B\ < 1 ) one 
should use Eq. (|33)l instead of ([TJ . Figure [3] shows the 
dependence g{M) given by Eq. (13"3"|) . The conductance 
varies smoothly from 2e 2 /h in the inverted regime (M < 
0) to zero in the ordinary insulator state with M > 0. 
Interestingly, there is a local minimum at M = 0. To 
demonstrate this analytically, let us examine Eq. (|33j) 
for M <C 7 and x = B M /\B\ <C 1. Using asymptotics 
Ti(x) ps Tr 1 / 2 ^ and r 2 (x) ps tt -1 / 2 , we expand: g(M) pa 
g(0) + Ag(M), where the correction Ag(M) is given by 



Ag(M) 

e 2 /^ 



7fl7L B M 
7 2 |B| 



7-R7L cM 2 



7 



2 4hv 2 \eB\ 



(35) 



It is quadratic in M reflecting the analyticity of the con- 
ductance in the limit M — > at finite B and 7. Therefore, 
on both sides of the transition the conductance initially 
increases with M, However, for M > the quadratic 
growth is followed by the suppression when the ordinary 
insulator regime with the gap M > 7 is reached. 



IV. CONCLUSIONS 

We have proposed a model for resonant magneto- 
transport via localized helical edge states in a quantum 
spin Hall system. To calculate the magnetoconductance 
we used the Green's function formalism allowing us to 
treat the topological and ordinary insulator regimes on 
equal footing. In the topologically nontrivial state, the 
magnetic-field dependence of the conductance has a res- 
onant (Lorenzian) line-shape explicitly depending on the 
orbital magnetic moment of the helical edge states. One 
could use this to extract the orbital g-factor of the heli- 
cal edge states in a given material. The necessary con- 
dition to test our predictions is a sufficiently large edge 
level spacing hv/L (see, Eq. (fMf ). For HgTe wells with 
v = 5.5 x 10 5 ms _1 and L = 1 /im the level spacing is 
about 2.5 meV (see, also, Fig. implying working tem- 
peratures up to a few kelvin. Note that with increasing 
B one should expect oscillations of the tunneling con- 
ductance that follow the zero-mode (B — 0) resonance 
because the states with nonzero momentum [n ^ in 
Eq. (|23[) ] will periodically cross the zero energy. Similar 
oscillations are expected in the gate voltage dependence 
of the conductance, which can be used to adjust the po- 
sition of the Fermi level in the band gap. 
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Appendix A: Electron confinement with symmetric 
boundary conditions 

We consider a system of finite length L (in the x di- 
rection in Fig. QJ,) subject to a perpendicular magnetic 
field B ([5]) and described by the equations 

ihd t ip+(r,t) = (Al) 
= (ihv[a x (d x + ieBy/ch) + a y d y ] - Ma z )ip + (r, t), 
ihBtil>-(T,t) = (A2) 
= {-ihv[a x (d x + ieBy/ch) + <J y d y \ + Ma z )ip-(r, t), 

for the operators ip± (r, t) of the Kramers-partner states 
of Hamiltonian At the ends of the system we impose 
the symmetric boundary conditions: 



i/>±(-L/2,y)=tl>±(L/2,y). 



(A3) 



Our goal is to show that the symmetric boundary con- 
ditions (IA3|) yield vanishing current density, 



j x (±L/2,y) = 0, 



(A4) 



and in this sense the particle is confined in the region 
|a;| < L/2. This becomes possible due to the specific 
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structure of the Dirac particle current density, 

3x(x,y) = ij}\(x,y)a x ip + (x,y) - ip]_{x,y)a x ip-(x,y)lAS) 

and the parity symmetry of Eqs. (|A1[) and (| A2|) . 

fli-(x,y;B)=C<r x il> + (-x,y,-B), \C\ 2 = 1. (A6) 

We now use this relation to transform the second term 
in Eq. (|A5[) as follows 

3x(x,y) = ipl(x,y;B)a x ip+(x,y;B) (A7) 
- ipl(-x,y;-B)a x il; + (-x,y;-B). 

Setting x — L/2 (or x = —L/2) and using the boundary 
conditions (| A3|) we have 

M±L/2,y) = ip* + (L/2, y; B)a x tp + (L/2, y; B) (A8) 
- ip\_(L/2, y; —B)a x ip + (L/2, y; —B). 



For zero magnetic field, B = 0, the two terms in Eq. (|A8I) . 
i.e. the currents of the Kramers partners, cancel each 
other out, yielding vanishing current density (|A4[) at the 
ends of the system. Moreover, the current density (IA8|) 
should also vanish in a finite magnetic field, B ^ 0, pro- 
vided that the system is in the quantum spin Hall regime. 
Indeed, the current density (|A8[) is odd in the magnetic 
field B and, therefore, could only characterize the Hall 
response. However, in the quantum spin Hall regime con- 
sidered in our paper a Kramers pair of counterpropagat- 
ing edge states cannot generate the Hall response. 
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